Significant Gamma Lines from Inert Higgs Dark Matter 
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One way to unambiguously confirm the existence of particle dark matter and determine its mass 
would be to detect its annihilation into monochromatic gamma-rays in upcoming telescopes. One 
of the most minimal models for dark matter is the inert doublet model, obtained by adding another 
Higgs doublet with no direct coupling to fermions. For a mass between 40 and 80 GeV, the lightest 
of the new inert Higgs particles can give the correct cosmic abundance of cold dark matter in 
agreement with current observations. We show that for this scalar dark matter candidate, the 
annihilation signal of monochromatic 77 and Z7 final states would be exceptionally strong. The 
energy range and rates for these gamma-ray line signals make them ideal to search for with the soon 
upcoming GLAST satellite. 

PACS numbers: 95.35.+d, 14.80.Cp, 98.70.Rz 



Introduction. — Recent measurements have established 
that nonbaryonic cold dark matter (CDM) makes up 
about one quarter of the Universe's total energy budget 
1 1 . Although a variety of candidates have been proposed 
2j , the nature of this dark matter still remains a mystery. 

One particularly attractive class of candidates is that 
of weakly interacting massive particles (WIMPs), as 
gauge couplings and masses of the electroweak symme- 
try breaking scale (i.e. presumably around 100 GeV) give 
the right order of magnitude for their thermal relic abun- 
dance to explain the dark matter. Moreover, electroweak 
precision tests indicate that the particle standard model 
Higgs boson should be light: run < 144 GeV at 95% con- 
fidence level, with a central value around 76 GeV 
However, in the standard model, the Higgs boson mass 
acquires quadratic ultraviolet divergences, requiring fine- 
tuning to keep the Higgs mass low if no new divergence- 
cancelling physics appears before the TeV energy scale 
(the hierarchy problem). Therefore, new physics can be 
expected to be found in upcoming high-energy experi- 
ments, such as CERN's Large Hadron Collider (LHC). 

The most studied scenario which provides both a sym- 
metry to cancel the quadratic divergences in the Higgs 
sector and a WIMP dark matter candidate is supersym- 
metry. In the minimal supersymmetric standard model, 
the lightest Higgsparticle is constrained to be lighter 
than ~135 GeV |4j, and although not excluded, some 
amount of fine-tuning [j| is actually needed to fulfill the 
experimental lower bound of 114 GeV This has mo- 
tivated several studies on how to theoretically allow for 
an increase of the Higgs mass both within supersymmc- 
try and other extensions of the standard model (see, e.g., 
[5j, [6j and references therein) . 

The Inert Doublet Model (IDM) @ considered in this 
Letter is a very minimal extension of the standard model 
- an added second Higgs doublet H 2 , with an imposed 
unbroken discrete Z2 symmetry that forbids its direct 
coupling to fermions (i.e. H2 is inert). In the IDM the 
standard model Higgs mass can be as high as about 500 
GeV and still fulfill experimental precision tests Q . Fur- 



thermore, conservation of the Z2 parity implies that the 
lightest inert Higgs particle (7J°) is stable and hence a 
good dark matter candidate @, @]. Although the IDM 
does not solve the hierarchy problem, but potentially 
only pushes the need for divergence cancelling physics 
beyond the reach of upcoming accelerator searches such 
as the LHC [1, H|, it has the advantage of providing a 
scalar WIMP dark matter candidate in a very minimal- 
istic way. An electrically neutral H° could therefore be 
used to represent a larger class of scalar dark matter, 
similar to how the neutralino in minimal supersymme- 
try often works as an archetype for supersymmetric dark 
matter. Actually, the origin of the IDM goes back to 
the 1970s [sl and has recently received new interest not 
only because of its potential to allow for a high Higgs 
mass and a dark matter candidate, but also forgener- 
ating light neutrinos and leptogenesis (see, e.g., [1 01 ] and 
references therein). This gives strong reasons to study 
detection prospects for the inert Higgs. 

A study of the relic density showed that H° can con- 
stitute all dark matter if its mass is roughly 10 — 80 GeV 



(or above 500 GeV if parameters are tuned) [111]. Di- 
rect detection studies show that the sensitivity of current 
experiments is far too low, whereas future experiments 
could cover more of the IDM parameter space [11]. The 
model can also produce observable signals at the LHC 
As shown in [ll|, indirect detection of the contin- 
uous gamma-ray spectrum might be reachable with the 
upcoming GLAST experiment [13]. However, this study 
was made for standard model Higgs masses of 120 and 
200 GeV which, although giving higher gamma rates, de- 
viates from the motivation for the model of a raised Higgs 
mass Q and lacks distinctive detectable features. 

An interesting aspect of the IDM is that the H° mass 
generically has to be below the charged gauge boson 
mass, since the relatively strong coupling to W + W~ 
would otherwise give a too low relic density to explain the 
dark matter. Virtual W bosons close to threshold could, 
on the other hand, significantly enhance loop processes 
producing monochromatic photons. We show here that 
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this is indeed correct and study the dark matter "smok- 
ing gun" line signals from the final states 77 and, when 
kinematically allowed, Z7. This, in combination with 
small tree-level annihilation rates into fermions, makes 
the gamma lines the most promising indirect detection 
signal. 

Calculations of monochromatic gamma lines for both 
spin- 1/2 [IH and spin-1 [3| dark matter annihilation 
have been performed earlier. The spin-0 dark matter can- 
didate discussed here has an even more promising gamma 
line rate. 

The Inert Higgs Model and its Constraints. — The IDM 
framework is an extension of the particle standard model 
with one additional Higgs doublet H2 and an unbroken 
Z2 symmetry under which Hi is odd, H2 — » —H2, while 
all other fields are unchanged (even). The potential for 
such a model is [(| 

V = m?|#i| 2 + V 2 2\H 2 \ 2 + Ailffil 4 + \2\H2\ 4 

+ XslH^l^l 2 + \i\HlH 2 \ 2 + \ 5 Re[(HlH 2 ) 2 ] , (1) 

where Hi is the standard model Higgs doublet, and \Xi 
and Xi are real parameters. 

Besides the standard model Higgs particle (h), the 
physical states derived from the inert doublet H2 are two 
charged states (H^) and two neutral: one CP-even (H°) 
and one CP-odd (A ), either of which is the dark matter 
candidate. For definitcness we have chosen H° as the 
lightest inert particle, although the role of H° and A 
can in general be interchanged. The masses are given by 



m 2 H± 



(4 + (A 3 + A 4 + \ 5 )v 2 
l4 + (A 3 + A 4 - X 5 )v 2 

f4 + x 3 v 2 , 



(2) 



where v = to^/V^Ai" is the vacuum expectation value 
for Hi. As usual, the W and Z masses determine v ~ 
175 GeV. No vacuum expectation value is allowed for H2 
because of the imposed unbroken Z2 parity. 

We apply the constraints in Q concerning vacuum sta- 
bility, pcrturbativity, precision tests, accelerator searches 
and naturalness [2J] while using direct detection and 
charged scalar mass bounds from ll|. We also impose 
m,H° + > mz to not be in conflict with data on the 
width of the Z boson. Furthermore, in agreement with 
observations, we constrain the H° abundance to 0.094 
< ftcDAih 2 < 0.129 [l|. The H° relic density calculations 
have been performed by interfacing the standard numer- 
ical packages FormCalc [l]| and DarkSUSY [lit . 

Gamma-ray Lines from Inert Higgs Annihilations. — 
For H masses below mw, only annihilations into 
fermions lighter than mjjo are accessible at tree level. 
The annihilation rate is given by 



N c na 2 m 2 



(i-^) 3/2 K. 
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where N c is a color factor (which equals 1 for leptons 
and 3 for quarks), yfs is the center of mass energy, a the 
fine-structure constant, mw the W boson mass, 0w the 
weak mixing angle, Th the decay width of h, and m/ the 
final state fermion mass. 

The heaviest kinematically allowed fermion state will 
dominate the tree-level annihilation channels. The con- 
tributions to the continuum gamma-ray spectrum from 
annihilating H° with masses below mw predominantly 
come from the secondary gamma rays produced in the 
fragmentation of the bottom quarks, mainly by the se- 
quential production and decay of neutral pions. Because 
of the harder gamma spectrum from the decay of r- 
lcptons they contribute significantly at the highest en- 
ergies, despite their much lower branching ratio. We use 



17J to calculate the shape of the continuum spectra. 

We now turn to the main issue in this Letter - the 
very important line signals from direct annihilation of H° 
pairs into 77 and Zj. These spectral lines would show up 
as characteristic dark matter fingerprints at the energies 



m H <> and m#o 



? /4m^o, respectively. The Z7 line 



might not be striclty monochromatic due to the Brcit- 
Wigncr width of the Z mass, but can still be strongly 
peaked. When the branching ratio into Z7 becomes sig- 
nificant the subsequent decay of the Z boson significantly 
contributes to the continuum gamma-ray spectrum. The 
full one-loop Feynman amplitudes were calculated by 
implementing the IDM into the numerical FormCalc 
package [H||. 

Four IDM benchmark models are defined in Table U 
where the two models III and IV are similar to those 
investigated in llj. Annihilation rates, branching ratios 



and relic densities for these models are given in Table HT1 
As an illustrative example, Fig. Q] shows the predicted 
gamma spectrum for model I. 

The spectral shape with its characteristic peaks in the 
hitherto unexplored energy range between 30 and 100 
GeV is ideal to search for with the GLAST experiment 
[lit . I n Fig- El this is illustrated by showing the predicted 
fluxes from a Ail = 10~ 3 sr region around the direction of 
the galactic center together with existing observations in 
the same sky direction. For simplicity, we assume a stan- 
dard Navarro, Frenk and White (NFW) density profile 
[lit for the dark matter halo in our galaxy (J x Afl <~ 1 
for AQ = 10~ 3 sr with the notation of [20j). Processes 
such as adiabatic compression could enhance the dark 
matter density significantly near the galactic center (see, 
e -g-j [H|), and we therefore allow our predicted flux to 
be scaled by a "boost factor" . 

The Energetic Gamma-Ray Experiment Telescope 
(EGRET) data, taken from [20j, set an upper limit for 
the continuum part of our spectrum. For example, for 
benchmark model II one finds that an optimistic, but 

boost of 10 4 might be al- 
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not necessarily unrealistic 
lowed. In such cases, there would be a truly spectacular 
77 line signal waiting for GLAST. However, to enable 
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FIG. 1: The total differential photon distribution from anni- 
hilations of an inert Higgs dark matter particle (solid line). 
Shown separately are the contributions from H H — > bb 
(dashed line), t + t~ (dash-dotted line) and Z'y (dotted line). 
This is for the benchmark model I in Tabled 



detection, boost factors of such magnitudes are not nec- 
essary. For H° masses closer to the W threshold the 77 
annihilation rates become even higher and in addition 
production becomes important. In fact, these signals 
would potentially be visible even without any boost at all 
(especially if the background is low, as might be the case 
if the EGRET signal is an galactic off-center source as 
indicated in (22)]). Also shown in Fig. [2] is the data from 
the currently operating air Chercnkov telescope HESS 
[23! ]. One may notice that future air Cherenkov tele- 
scopes with lower energy thresholds will cover all of the 
interesting region for this dark matter candidate. 

Finally, we have made a systematic parameter scan 
for rrih = 500 GeV, calculating the cross section into 
gamma lines. The previously mentioned constraints al- 
low us to scan the full parameter space for dark matter 
masses below the W threshold of 80 GeV. The depen- 
dence on m#± and A2 is small, and we set these equal to 
mjjo + 120 GeV (to fulfill precision tests) and 0.1, respec- 
tively. Importantly, one notes that the right relic density 
is obtained with a significant amount of early Universe 
coannihilations with the inert A particle. The resulting 
annihilation rates into 77 and Z7 are shown in Fig. [3] 
The lower and upper m H o mass bounds come from the 
accelerator constraints and the effect on the relic density 
by the opening of the W + W~ annihilation channel, re- 
spectively For comparison, we show in the same figure 



^ -4 
> 

CP 

O -5 

I 



-6 - 



3L - 7 

® -8 

(m ^■ 

3 -9 
o 

J -10 



IDM: NFW, An~10" d , cr E =7% ' 


t 


np\/ hnnct~m 4 
OU OfcJ V , UUUbl ~ I u 


EGRET:An=2x1(T 3 








■ y< \ 

, &>«i / 


.70 GeV, boost -100 : 


X 






12 3 4 

Log(£ 7 [GeV]) 

FIG. 2: Predicted gamma-ray spectra from the inert Higgs 
benchmark models I and II as seen by GLAST (solid lines). 
The predicted gamma flux is from a AO, = 10 -3 sr region 
around the direction of the galactic center assuming an NFW 
halo profile (with boost factors as indicated in the figure) and 
convolved with a 7 % Gaussian energy resolution. The boxes 
show EGRET data (which set an upper limit for the contin- 
uum signal) and the thick line HESS data in the same sky di- 
rection. The GLAST sensitivity (dotted line) is here defined 
as 10 detected events within an effective exposure of 1 m 2 yr 
within a relative energy range of ±7%. 



the corresponding annihilation rates for the neutralino 
(x) within the minimal supersymmetric standard model. 
The stronger line signal and smaller spread in the pre- 
dicted IDM flux are caused by the allowed unsuppressed 
coupling to W pairs that appear in contributing Feynman 
loop diagrams. 

Summary and Conclusions. — In this Letter, we have 
investigated the gamma-ray spectrum from the annihi- 
lation of the inert Higgs dark matter candidate H°. In 
particular, we have focused on its striking gamma lines 
which arise at the one-loop level and produce an excep- 
tionally clear dark matter signal. 

The gamma line signals are particularly strong for this 
scalar dark matter model mainly for two reasons: (1) The 
dark matter mass is just below the kinematic threshold 
for W production in the zero velocity limit. (2) The 
dark matter candidate almost decouples from fermions 
(i.e., couples only via standard model Higgs exchange), 
while still having ordinary gauge couplings to the gauge 
bosons. In fact, these two properties could define a more 



TABLE II: IDM benchmark model results. 
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IDM benchmark models. 
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FIG. 3: Annihilation rates into gamma-ray lines 2ucr 77 (up- 
per band) and voz-y (middle band) from the scan over the 
IDM parameter space. For comparison the lower-right region 
indicates the corresponding results within the minimal super- 
symmetric standard model as obtained with the DarkSUSY 
package Xl d /. 



general class of models for which the IDM is an attractive 
archetype. Despite small H° annihilation cross sections, 
coannihilations with A in the early Universe can bring 
the relic density of H° into the correct range. The com- 
bination of the low (or even completely vanishing) tree- 
level annihilation rates today and the strong loop-level 
processes, due to vmsuppressed couplings to virtual W 
bosons close to threshold, make the gamma lines the far 
most dominant feature in the resulting gamma spectrum. 

Absolute gamma-ray fluxes arc, unfortunately, still 
hard to predict due to the uncertainties in the structure 
of dark matter halos. It might eventually be the spec- 
tral shape that enables a separation of a dark matter 
signal and the background, in which case a gamma-line 
would be a striking feature. We have shown that such sig- 
nals in the IDM are promising features to search for with 
the GLAST satellite and with future air Cherenkov tele- 
scopes. One should bear in mind that the best prospects 
for detection might not be in the direction of the galactic 
center, but rather for other sources, such as dwarf galax- 
ies, smaller dark matter clumps or in the extragalactic 
gamma-ray radiation, where the background is lower. 

M.G. thanks Jan Conrad and Alexander Sellcrholm for 
useful discussions on the GLAST performance. M.G. and 
E.L. thank Thomas Hahn for technical assistance with 
FormCalc. L.B. and J.E. are grateful to the Swedish 
Science Research Council (VR) for support. 



t 
t 
§ 

[1 
[2 
[3 
[4 
[5 
[6 
[7] 
[«] 
[9 
[10 

[11 

[12 
[13 

[14 
[15 

[1(3 

fi- 
ll* 
[19 
[20 
[21 
[22 
[23 

[24 



Electronic address: 
Electronic address: 
Electronic address: 
Electronic address: 
D. N. Spergel et al 



michael@physto.se 

erik@physto.se 

lbe@physto.se 

edsjo@physto.se 

, Astrophys. J. Suppl. 



Ser. 170, 377 



Silk, Phys. Rept. 405, 279 



Espinosa and I. Hidalgo, Nucl. Phys. 
[hep-ph/0607279] . 

Rev. D 18, 2574 



(2007) [astro-ph/ 0603449 . 
G. Bertone, D. Hooper and J. 
(2005) [hep-ph/0404T75] . 
LEP Electroweak Working Group, status of March 2007, 
|http://lepewwg.web.cern.ch/LEPEWWG/| 
W. M. Y ao et al. (Particle D ata Group), J. Phys. G 33, 
1 (2006) ( [http://pdg.rbl.govl 

R. Barbieri, L. J. Hall, Y. Nomura and V. S. Rychkov, 
Phys. Rev. D 75, 035007 (2007) [hep-ph/0607332] . 
R. Barbieri, L. J. Hall and V. S. Ry chkov, Phys. Rev. D 
74, 015007 (2006) [hep-ph/0603"l88] . 
E. Ma, Phys . Rev. D 73, 077301 (2006) 
[hep-ph/060T225] . 
J. A. Casas, J. R 
B 777, 226 (2007) 

N. G. Deshpande and E. Ma, Phys 
(1978). 

E. Ma and U. S arkar, Phys. Lett. B 653, 288 (2007) 
|arXiv:0705.0"074l [hep-ph]] . ' 
L. Lopez Honorez, E. Nezri 
M. H. G. T ytgat, JCAP 
[hep-ph/06i2275] . 

N. Gehrels and P. Michelson [GLAST Collaboration 
tropart. Phys. 11, 277 (1999). 

L. Bergstrom and P. Ullio, Nucl. Phys. B 504, 27 (1997) 
[hep-ph/9706232l ; P. Ullio and L. Ber gstrom, Phys. 
Rev. D 57, 1962 (1998) [hep-ph/970"7333] ; L. Bergstrom, 
P. Ullio and J. H. B uckley, Astropart. Phys. 9, 137 (1998) 
[astro-ph/9712318] . 

L. Bergstrom, T. Bringmann, M. Eriksson & M. Gustafs- 
son, JCAP 0504, 004 (2005) |hep-ph/041200Tj . 
T. Hahn and M. Perez- Victoria, Comput. Phys. Com- 
118, 153 (1999) hep-ph/9807565 ; T. Hahn, Com- 



J. L. 
0702. 



Oliver and 
028 (2007) 



As- 



mun. 118, 153 (1999) [hep-ph/980756 
put. Phys. Commun. 140, 418 (2001) 
P. Gondolo, J. Edsjo, P. Ullio 
M. Schelke and E. A. Baltz, JCAP 
[astro-ph/0406204| . 

T. Sjostrand, S. Mrenna and P. Skands 
(2006) [hep-ph/0603T75] . 
G. Zaharijas and D. Hooper, Phys. Rev 
(2006) 



hep-ph/00 12260] , 
L. Bergstrom, 
0407, 008 (2004) 

JHEP 0605, 026 



D 73. 103501 



astro-ph/0305075 



astro-ph/0603540 . 
J. F. Navarro, C. S. Fre nk and S. D. M. W hite, Astro 
phys. J. 490, 493 (1997) [astro-ph/9611107l . 
A. Cesarini, F. Fucito, A. Lionetto, A. Morselli and P. Ul 
lio, Astropart. Phys. 21, 267 (2004) 
G. Bertone and D. Merritt, Phys. 

(2005) [astro-ph/056T555] . 

D. Hooper and B. L. Dingus, Phys. Rev. D 70 
(2004) [astro-ph/0210617| . 

F. Aharonian et al. [H.E.S.S. Collaboration] 
Rev. Le tt. 97, 221102 (200 6) [Erratum-ibid. 97, 

(2006) ] [astro-ph/0610509] . 

Because of arbitrariness in defining naturalness, we relax 
this constraint by a factor of 2 compared to [(J . 



Rev. D 72, 103502 



113007 



Phys. 
249901 



